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Meta Heuristic methods have made a deep impact in the area of optimization 
in different streams of engineering. The performance of these algorithms is 
of importance because the hardware implementation of these algorithms is to 
be carried out for different engineering applications. As an important 
application in High Voltage DC (HVDC) transmission and Industrial Drives 
the multilevel inverter fault diagnosis is carried out using the different meta¬ 
heuristic methods with Neural Network as the decision making algorithm. 
The optimization of the weight and the bias values in the neural network 
diagnosis system is carried out in order to analyze the performance by means 
of comparing the Mean Square Error (MSE) while the Neural Network is 
getting trained for different fault conditions in the multilevel inverter. Matlab 
based implementation is carried out and the results are tabulated and inferred 
for a Multilevel Inverter fed from the Photovoltaic power generation system. 
In order to increase the robustness of the fault detection, with renewable 
energy based power generation system as the source for the Multilevel 
Inverter, the feature extracted from the multilevel inverter are positive, 
negative and zero sequence voltage along with the THD of the output 
voltage. The optimization algorithm used is Particle Swarm Optimization 
(PSO), Cuckoo Search Algorithm(CSA), Genetic Algorithm(GA) and Tabu 
Search Algorithm (TSA) 
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I. INTRODUCTION 

The The usage of multilevel inverter is expanding day by day as the HVDC and the Industrial 
Drives need performance enhancement in the quality of power that has to be delivered to the industry and the 
transmission systems.Cascaded H bridge being the advanced multilevel inverter topology,many literatures 
are using the Cascaded topology. Optimized Artificial Neural Network based implementation for fault 
detection of multilevel inverter is carried out in [1-3].Investigation of different fault mode behavior is carried 
out for different types of faults [4]. The park’s transform vector pattern is used for the fault diagnosis in the 
voltage source inverter fed AC drives [5] [6]. The classification of the fault wasdependent on the load 
characteristics as current is considered for fault detection.The different fault detection methods are reviewed 
[7], voltage histogram based fault detection is carried out in [8-9].The wavelet analysis of the output voltages 
and input DC currents for fault diagnosis is carried out in [10-11].Since the current waveform would be a 
problem when large variation in the loads would lead to misinterpretation the wavelet analysis based fault 
detection is carried out [12]. Output voltage frequency analysis [13], artificial intelligence based fault 
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detection and reconfiguration [14-15] are also proposed which are complex to implement. The different 
meta-heuristic based fault identification and the remedial action are also carried out in literature [17-21]. A 
time analysis of the zero current occurring while the zero voltage is occurring in a multilevel inverter, which 
is controlled by the rotating level shifted pulse width modulation is carried out [22]. Probabilistic Neural 
Network for learning the fault detection using the wavelet as the feature extraction is carried out in [23].This 
paper attempts to carry out a performance analysis of the different meta- heuristic method on the diagnosis of 
the open switch fault in the cascaded multilevel inverter. The fault conditions are trained on the ANN with 
the radial basis function as the learning function. The parameters like the MSE is optimized by the use of the 
Meta heuristic methods like PSO,GA, TSA and CSA.A robust algorithm for fault detection is a need, as in 
the previous literature; it is obvious that the fault prediction has created more misinterpretation. Thus the 
learning based algorithms have to be implemented which would facilitate the robust fault detection. Thus the 
radial basis function based ANN implementation is chose. The training of the RBE based ANN is carried out 
by optimizing the training parameters comprising both weight and the bias values in the neural network 
learning process. This paper which is the improvement of the previous works [1-3] is organized as follows, 
the Section -II would have the description about the training of the ANN with the radial basis function. 
Section -III discusses about the optimization algorithm used in the ANN parameter estimation, Section- 
IVexamines the results and discussion with inference of the different meta heuristic methods, and it is 
followed by the conclusion and references. 


2. FORMULATION OF PARAMETER ESTIMATION USING META HEURISTIC METHODS 

The formulation of any meta-heuristic method would include the identification of parameter or the 
independent variables, initial number of population of the variables chosen, objective function, vector 
function of the meta-heuristic method. The parameters already chosen in the previous section are the weight 
and the bias value of the radial basis function regression. The Figure 1 depicts the overall implementation of 
the parameter estimation algorithm that we implement in this paper. 



Figure 1. Overall parameter estimation block diagram 


The topology of the three phase Cascaded multilevel Inverter is as shown in Figure 2. The objective 
function is the reduction of the MSE and RMSE during the training of the RBFBANN.The velocity update 
function of the PSO is 
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V;(? -h 1) = wv,.(0 -h c^r^{Xi{t) - x,.(0] + c^r^{g{t) - x,.(0] 

A 

where V, (t + 1) is the new value of the weight or the bias value , Xi{t) is the best weight or bias value in 
the current iteration, Cj and C 2 are constants with values near 2, g{t) is the weight or bias value which is 
the global best in all the iteration calculated so far, ^ is the inertial coefficient which is a constant. 



Figure 2. Seven level three phase cascaded multilevel inverter topology 


The CSA algorithm is inspired by the way the Cuckoo would act as an intruder to the crow’s nest to 
lay eggs in the erow’s nest. The best nest would allow the egg to grow but the worst nest would throw away 
the eggs. The velocity update function for the CSA algorithm is as given [18], The probability of the host 
bird to find the intruder egg is Pa. And the new nest can be found by the following formula, 

X.(t + l) = Xi(t) + alevy(/i) (2) 

where Xj (f -I- 1) is the possible solutions , X.{t) is the previous solution, is the step size which 

controls the amount of search scale, is the number of nest that has to be abandoned in the current iteration. 
The Cuckoo would lay the eggs by a randomness of the Levy Flight, which is defined by the 
following function. 


levy(X) = 


r(l-H A) X sin(nA / 2) 
r((i -H ;i) / 2) X A X 


where is a constant (1< 3) and is the gamma function which is defined by the integral value , 


r(x) = Je'rtf 

0 
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as the velocity update function of GA the following process would be calculated.The parent selection for the 
population generation is carried using the Roulette Selection method as defined below[3], 

S, =—^— 

Sum(Y) 

S1 = 0 

J =1 


Where Y is the set of all fitness values, is the ith fitness value. Let SI be some random value and S2 and 
J are variables. The Roulette selection algorithm can be defined as follows, 

While -S' 2 > 51 the following steps has to follow 

Find^2 = S2+S(J) 

Increment J 
End 

Once the while loop is becoming false the J value is taken and the population at the J th position is 
selected as the parent. 

The arithmetic crossover is used to create the new offspring from the two parents selected using the 
Roulette selection. Two off springs are generated using the following formulae. 

Offspring l = a* parent 1 -i- (1 a)* parent 2 
Offspring 2 = (1 a) * parent l+a* parent 2 

Where, a is the random weighting factor. In order to maintain genetic diversity between the off springs the 
mutation operation is carried out[3]. For the calculation of the mutation the maximum and the minimum 
weight values are taken from the a values and it is mentioned in the below equations as ‘maxweight’ and 
‘minweight’. The mutation operation is defined in the following formulae , 

= 1 2*randnum 

Where , fctndnum js the random value. 

Using the above value the inverse error function is calculated and utilized in the following formulae, 

Z =inverseerrorfunction( )*2“^ (3) 

The offspring is calculated by using 

Offspring = Z * + parent 


{msixweight mmweight ) 
10 


reduces its value for every iteration ' 

The new populated value from the genetic algorithm is used to be applied in the objective function. The 
updation of the parameter in the TSA would occur by finding the neighbors in the solution space. TSA 
develops the new solution space using the memory structures using the equation (10) 

N _list(i,:) = sbest * rand* S(randi(pop _size,l),:) (4) 

The solution space for the new iteration is developed by using the above equation where 
N _ list {ip) is all the neighbors that would be used in the next iteration wAsbest is the best solution that 

occurred in the previous iteration and S is the set of all the neighbors which is the memory. This list of new 
neighbors is called the Tabulist, which is used for finding the solution by evaluating the objective 
function.The comparison of all the above meta-heuristics based implementation is carried out and the results 
are tabulated. The comparative analysis of the different methods are analyzed and inferred..By different 
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experiments it is come to a conclusion that RBF is a universal approximating function and classifier, thus it 
happens to be used in this paper. 

3. RESULTS AND ANALYSIS 

The Matlab based implementation is carried out on the 7 level Multilevel inverter and the execution 
time of the RBFBANN without the optimization, with PSO,GA,TSA and CSA is carried out and the results 
are tabulated and inferred for the performance analysis of the different methods.The results that are obtained 
from the basic RBFBANN are as given below for the same kind of graphs would be obtained for all the 
methods thus not repeated. The Figure 3 depicts the positive voltage sequence that is given as the 
instantaneous input to the input nodes of the neural network. Similarly Figure 4,5 are the Negative and the 
zero sequence voltages of the inverter respectively. The Figure 6,7 and 8 are the positive, negative and zero 
sequence phase respectively of the inverter that is fed to the input nodes of the Neural Network. For the time 
secs 0 to 0.25secs and 0.75 to Isecs there is a phase difference in angle calculation. So these three is taken as 
another fault classifier input for neural network. The difference is zoomed and shown in the above figure. 
These are calculated from equation from the results it can be seen that the basic radial function is used to 
reduce the mean square error (10-29). Therefore RBF is used in neural networks and trained with PSO to 
reduce the average square error. And that gives a very less error value, around zero. 

3re has real and imaginary value real value is magnitude and imaginary is angle. 

Figure 9 is the important graph, which depicts the instantaneous THD value of the inverter at 
different switch failure. The Figure 9 shows the voltage THD values of Ql, Q2, Q3 and Q4 faulted at each 
0.25 sec respectively. So this THD value is given as input and target is used as the switch number 1 to 4 in 
each 0.25 sec. After learning using this faults can be identified and remedial actions can be taken. The 
comparative analysis of the execution time of the different algorithms with different iteration number is given 
in Table 1. By using cuckoo search algorithm, the error minimizes to zero consequently more accurate than 
PSO. So RBF neural network trained with cuckoo search algorithm is performing better comparatively. 
According to the MSB these algorithms are getting the lowest of all the algorithm is given by the cuckoo 
search algorithm. The Table 2 shows about the execution time of the Matlab code on the computer system 
with 8 GB RAM and 2.4Ghz Intel i5 processor. The execution time of the CSA algorithm is comparatively 
less and thus it has the advantage of being implemented on the hardware. Hardware feasibility is more for the 
CSA algorithm. 


= -* (Y^ + a * Vij + * 

1 

* (Va + a * V,, + * 

1 

'^2 (Va + * Vi, + a * 

= + Where, 
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Figure 3. Positive voltage sequence 
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Voltage at - seq 



Figure 4. Negative voltage sequence 
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Figure 5. Zero voltage sequence 



Figure 6. Positive sequence phase with and without fault 
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Figure 7. Negative sequence phase with and without fault 
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Figure 8. Zero sequence phase with and without fault 
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Figure 9. THD A-phase 
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Table 1. Mean Square Error Comparison of Neural Network trained with and RBE Neural Network trained 

with PSO and CSA 


Run no. 

RBF with PSO 

RBE with GA 

RBF with TSA 

RBF with CSA 

1 

2.11E-28 

2.90E-25 

6.44E-29 

6.64E-29 

2 

1.44E-28 

1.33E-25 

7.95E-29 

7.13E-29 

3 

1.03E-28 

1.23E-25 

1.87E-28 

1.47E-28 

4 

1.02E-28 

1.12E-25 

6.89E-29 

6.38E-29 

5 

1.21E-28 

l.OlE-25 

1.88E-28 

1.12E-28 

6 

2.11E-28 

l.OOE-25 

1.88E-28 

1.13E-28 

7 

7.12E-29 

3.12E-26 

8.65E-29 

8.74E-29 

8 

1.84E-28 

1.94E-26 

7.16E-29 

7.11E-29 

9 

1.04E-28 

1.74E-26 

7.635E-29 

7.11E-29 

10 

9.04E-28 

5.04E-26 

8.73763E-29 

8.36E-29 


4. CONCLUSION 

A highly robust and quick fault diagnosis system in a three phase 7 level Multilevel inverter using 
the optimized RBEANN comprising meta-heuristic methods like PSO,GA,TSA and CSA. The single switch 
open fault for all the switches are trained using different algorithms and the results are discussed and 
tabulated. The fault identification after getting trained is tested in different condition and was found to be 
robust. According to the performance measure of MSE and RMSE the CSA algorithm performed better with 
the lowest MSE and RMSE and execution time was found less compared to other algorithms taken. Thus due 
to the reduced MSE and RMSE the CSA would be preferred for the hardware implementation. 
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